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FOREWORD 

'Olio  publication  vac  prepared  under  contract 
by  the  UNITED  STATES  JOINT  PUBLICATIONS  RE¬ 
SEARCH  SERVICE,  a  federal  government  organi¬ 
zation  established  to  service  the  translation 
and  research  needs  of  the  various  government 
departments . 
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'  L!  :;FOR  AUTOMATIC  COMPENSATION  OP  ACCELERATIONS  DUE  TO  GRAVITY  j 


‘  7  /  [Following  i6  a  translation  of  an  article  by  V.  A.  Bodner  uif.cf 

-  V.‘  P.  Seleznev  in  the  Russian-language  periodical  ■  .d*:; 

31  fevastiya  Akaaaaii  NhUk  SSSR , 3 OTB .  Energetika  i  Avtonatika, 

'  o :  '■  ^Bulletin  of  the  Academy  of  Sciences  USSR,  Department  •  of  .x  i 
i  y  Technical  Sciences/  Power  Engineering  and  Automation),  No.  1, 

1  Moscow,  January- February  I960,  pages'  76-85 .  o=sb  x^cX-.r^-r 

■  .  •  ‘  .  .  ■  •  iM-  f  *- '  ■  - '  oi>.  j'X  „Cj*  /.*••  #.*•;./■>*  0  !>;' j  C:> 

o  '•/Controlled  motion  of  a  body  Is  possible  in  base  the  coordinates' 

’•  fcf  its  '  position,  velbcity,  'ahd  6ther  navigational  elenents  in  respect 
r'*  to  a  determined  coordinate  system  are  known.  The  continuous  determina¬ 
tion  of  navigational  elements  is  effected  by  navigational  systens  which 
r;f : be  based  on  the  use  of  the  'properties  of  physical  fields::  (the  gravi- 
^  tdtlOnal  field,  magnetic  arid  electric  'fields,  etc...  Or  by  bearings  on 
celestial  bodies. Nonperturbed  systens,  which  are  frequently  ca.ibir-.ed 
r/ith  astronomical  systens/  occupy  a  special  place  anong  these  systems. 

An  astroinertial  system  in  which  the  favorable  properties  of  an  astro- 
•'  nailcai  systen  add  a  noripefturbed  system  are  combined  is  examined  in 
!  '-this  article  [1]  . is  However/  the  author  has  dot  taken  into  account  the 
peculiarities  of  nonperturbed  systems  with  three  channels  of  automatic 
'-compensation  and  has  not  evaluated  the  errors  of  such  systems. 

■in:i;x..:u>o:  ijkis  article  contains  a  discussion'  of  certain  problems  of  the 


/.theory  of  hydroinert ial  systems  operating  in  three-dimensional  space 
•• '  ‘  Jand  an  analysis  of  netho&i cal  errors .  :.?/)  ■  x? 

1.  Principles  of  the  Construction  of  nonperturbed  Systems. ^ 

-  The ‘navigational  's'lenents  in  a  oionperturbed  system  are  determined  on 
the  basis  of  tnc  integration  of  accelerations  measured  by  three  ae- 
.  celeroieters  with  mutually  perpendicular  axes.  The  accelerations  are 
rX  measured  in  respect  to  an  inertial  coordinate  system  which  has  been 
selected  in  such  a  manner  as  to  have  no  angular  velocities  relative 
-rt0'  stellar  space  and  whose  ^origin  can  move  'under* the  action  of  gravi¬ 
tational  forces .  The  accelerometers  measure  only  the  accelerations 
caused  by  the  action  of  forces  of  drag  and  the  thrust  of  the  engines, 
but  not  of  gravitational  forces.  "  .1:1 


1 


Lst  us  select  a  coordinate  systen.  The  origin  of  the  coordi- 
riiates  tdf  the  inertial  systen  can  coincide  with  the  center  of  inertia 
0  -  bf  'hny  .celestial  body  which  is  moving  in  a  gravitational  field.  In 
flight  near  the  earth,  an  equatorial  system,  two  axes  of  which  coin¬ 
cide  with  the  plane  of  the  equator  and  the  third  axis  with  the  axis 
8‘i 9oif rbthtiOh  of  the  earth,  'can  be  utilized-. as,  the  inertial  -coordinate 
systen.  The  chief  shortcoming  of  such  a  systen  is  the  change  in  the 
position  of  .  its Jaxesi  relative  to  stellar  space  due  to  the  rotation 
of  the  earth. 

For  interplanetary  flights  it  is  convenient  to  make  use  of  an 
elliptical  systen  of  coordinates  Xq,  y0,  zQ  (Figure  l)  with  its  origin 
in  the  center  of  the  sun.  Two  axes  (xq,  Fq)  lie  in  the  plane  of  the 
earth's  orbit  (ecliptic),  and  the  third  axis  is  perpendicular  to  this 
plane.;  tOiie  Of. 'the /axes:,  (xq).  can  be  directed  toward  the -fir ?t -magni¬ 
tude  star  Spica,  which  is  located  near  the  plane  of  the  fcliptic.  Die 
position  of  the  flying  apparatus  in  such  a  coordinate  systen -is  deter- 
nined^by ^rectangular  coordinates  x,  y,  z.  It  is  also  possible  to  make 
use  of  a  spherical  coordinate  systen  with  the  sane  origin,  positions 
in  which  are  determined  by  the  astrpnonical  latitude  v.ganto'.ond-  lohgi> 
tude  Inrihrln.  with  radius  vector  R. 

!  'Jy\ The  ecliptic:  coordinate  system  possesses  great  stability .  Thus, 
-for  example',  the -change:  in  the  astronauical  latitude  amounts  to  0.07 
nicroradians  per  year,  and  the  change  in  the  longitude  is  not  nore  than 


i-’:0.1  microradians; per  year.  .1. •. •••  ‘j..  f„v>.L. 

>-:■>)  Ihe"  following  forces  are,  acting  on  $ie  flying  apparatus  :  thrust, 

•  ‘  drag,aahd  gravity.  -  The  vector  of  the  first  two  forces  is,  designated  by 
>  f  and  the  vector  of  the;  force  of  gravity  by  G.  The  components  of 
• '  these  forces  on  the  axes  of  the  rectangular  coordinate  system  XQy0Zo 
• -'are  equal  to  Fx,  Fv,  Fz:  and  ^  <  %  respectively  .  let  ust  use 

X,  ’y, '  z  and  x’ ,  y,'i>  z»  to  designate  the  coordinates,  and  the  velocities 
of  the  center  of  mass  -of,  the  flying  apparatus  in  the  selected,  inertial 
:  r  coordinate -system  XQyQZQ. M'"::.. 

•  Expressions  for  the  components  of  the  acceleration  of  the  cen¬ 
ter  of  mass  of  the  flying  apparatus  can  be  obtained  from  the  equations 

of  notion  '  •  .  :  : ■  .  !'-V\  ’’  . 

nx"  =  Fx  -*•  Gx,  my"  -  Fy  —  Cfy-,  ,  hz  s  Fz  —  Gj-  (1*1) 

where ‘  n  is  the  mass  of  the  flying  apparatus.  .  I 

;  ■ :  if  ve  make  use  'of  the  expressions  for  the  accelerations  of  the 

flying  apparatus 

■  %o  *  Si,  a.„  5  £y,  azo  =  ^  (2.^) 
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' 'and  introduce  the  corresponding  accelerations  from  the  gravitational 
forces  '  '•  '-t-"'  v*  i.:-\  <  o::  ■  ■...■■2  >; 

:  6*  =  •  •  '  «z =  S&>  -  il*3) 

n  n  .  > 

then  we  find  from  equation  (l.l), 

axo  =  x"  +  &x>  =  y"  * 


1  >. 


=  b"  +  8a 


(1.4) 


'  Accelerometers  whose  sates  of  sensitivity  coincide  with  the 
axes  3^y0z0  will  he  used  to  measure  the  accelerations  of  (1.4).  If 
notion  occurs  only  under  the  action  of  gravitational  forces  (F  »  0), 
then  the  accelerations  of  (1.4)  vanish  and  the  aoceleroneters  will 
not  show  anything;:  "':;'; ;:r  ^  -  -f'  ■  f-  •■■■■■  r- • 

A  hydroinertial  systen  should  include  a  stabilized  platform, 
three  aoceleroneters  with  mutually  perpendicular  axes  of  sensitivity, 
integrators,  a  couputer,  and  feedback  e’lenents  (Figure ;2).  0iThe  axes 
of  sensitivity  of  the  accelefoneteis  ford  ah  orthogonal  Coordinate 
systen  xyz  which  should  be  parallel  with  the  inertial  coordinate  sys¬ 
ten  •  Due  to  the  drifting  of  each  of  the  gyroscopes,  however, 

the  xyz  systen  of  coordinates  deviates  fron  the  XQygZQ  systen  by  the 
^  angles'  alpha,  beta,  gonna  (Figure  3)* ; '  i  '  f"’;;  ""r;  .  • 

1  '  ,u':  ^  The  integrators  are  "necessary 'for  integrating  the  accelerations 
and  for  obtaining  the  velodity  And  coordinates  of  position.1  The  con- 
"puter  and  .the  feedback  systen  produce  signals  to  compensate  '‘for  ac¬ 
celerations  fron  gravitational  forces.  "  A  shortcoming  of  the  hydro¬ 
inertial  systen  is  the  fact  that  it  cannot  independently  eliminate 
the  angles  alpha,  beta  and  ganda.  .External  information  is  needed  to 
elininate  then. 

S.  Equations  of  Motion  of  a  Nonperturbed  gysten.  In  carder  to 
derive  the  equations  of  notion  of  aft  inertial  systen  it  is  necessary  • 
to  mice  use  of  the  relationship  between  the  components  of  the  accel^ 
^rations  along  the  XqYqZq  and  the  xyz  axes.  The  components  of  the 
accelerations  ax,  ay,  az  measured  by  the  accelerometers  are  not 
equal  to  the  components  of  acceleration  of  the  tenter  of  the  mss  Of 
the  flying  apparatus  a^,  a^,  az0,  The  relationship  between  these 
components  is  given  in  the  table  of  direction  cosines: 

"■  .aU'X/ffV v'V; ’ ' 

x  ; .  ..cos alpha  cos  beta  ,  /' 

y  sin  alpha  sin  gamma  -  cos'  alpha  sin  beta  cos  gamma  .  ..  : 

z'.  sin  alpha  cos  gatsna  cosa&jiipha  sib  beta  ..sin  ganm  ;  '-labs  beta  sin 

si;  sin  beta 

y  VTf,,T  ...  cos  beta  cos  ganrn  .  Mr  '  i'y 

z:;. ^  -  cos  beta. sin. ganda  ;;;i ; ; ;s;  .• 

i*  ;V;'‘  ,yr.  CT  ^  ^  ,  ... 

x  :  •  -  sin  alpha  dos  'beta  '/  "  *  ■  yg*‘  •’  '  "s.i- 

y,  cos  alpha  sin  gamma  sin  alpha  sin  beta  sin  slpha  — . 
zUoSs  alpha  cos  gamma-  sin  alpha  sin  beta  sin"  gamma 

If  we  make  use  of  the  fact  that  under  real' conditions  the  angles 
alpha,  beta,  and  gamma  are  snail,  then  the  sines  can  be  replaced  by  the 
angles  thenselves,  the  cosines  can  be  assumed  to  be  equal  to  unity,,  and 
tfie'  products  of  the  angles  can  be  neglected ' as  quantities  ;of  the ‘second 
order  of  smallness .  Then  the  relation  between  the  components  of  ac- 


i.c;jceleration' along  the  xyz  and  the  XoyQZQcaxesvwIllbbe : • 

ax  8  ax0  axr 

SLy  si  ay©  “%•  Delta  Oyj. 

as  *  ez0  Del'ba  azr 
-  3  - 


(2.1) 


where 

Delta  exr  .«  OyQ  Beta  ***  aZQX,  Delta  ayj.  ■  aZQ  —  aXQ  Beta,  Delta 
&zr  B  £xOx  OyQ  Gorjpa  : .  v  (2.2) 

(  v  :  ')  x  :  “•  ..  ■  i.'.  /■  ..  '• 

.  r  v  The  errors  Delta  axg  Delta  ayg  and  Delta  axg  occur  due  to  ;j. 

drifting  of  the  gyroscopes,  as  a  resuli  of  which  the  axes  of  sensi¬ 
tivity  of  the  acceleroneters  do  not  coincide  with  the  axes  of  the 
inertial  coordinate  systen. 

",  ,®ie  output  signals  of  .the. acceleroneters,  as  nay  be  seen  from 
the  strucjiilral  diagram  in  Figure  4 ,  pass  through  the  integrators 
with  transfer  numbers  k^,  kn^,  .^iz>.'^2x>  ^2y>  k2a  into  the  conputef 
which, computes  the  compensating  signals  gxc,  Syc^rSzc^'t°  c<xqpensate 
for  accelerations  from  gravitational  forces.  :  ,  ,,  ,  . 

The  coordinates  of  position  Sy,  Sz  and  the  velocities  of 
itho  center  of  nass  S^,  Sy*,  Sz'  of  the  flying  apparatus  are  calculated 
las  a  result  of  the  operation  of  the  inertial  systen.  ■  ■, .... 

.....  Fron  the  structural  diagroa  (Figure  4)  we  obtain  the  following 
equations  of  notion  of  each  of  the  channels  of  the  inertial  systen; 

l(*z  A-  xo'i  ~  .4WS$*"  (2.3) 


+  V  j  ;  ■V  ¥o  =  ^y 


g  )^la  -5-  z 
^zr J -  * 

p  - 


B’l 


It 


2z  -V 
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where  xQ,  y0,  L  and  Xq’,  y0',  zQ'  are  the  initial  values  Of  the  co¬ 
ordinates  and  the  velocities  of  the  flying  apparatus. 

It  can  be  seen  fron  the  structural  diaagran  (Figure  4)  that 
the  signals  fron  the  output  of  the  first  integrators  are  the  compon¬ 
ents  of  the  absolute* velocities  of  flight,  while  the  signals  ar  the 
output  of  the  second  integrators  are  the  components  of  the  distance 
*fc  IT  01V6  HgcL  • 

Let  us  transform  the  equations  of  (2.3),  for  which  purpose  we 
shall  select  transfer  nunbers  of  the  integrators  fron  the  condition^ 

kix^  =  VC2y  S  *1**28  =  1*  We  Sha11  find  •: 


Delta  x"  =  Delta  g^.  -?  Delta  a^  (2.4) 

Delta  y"  -  Delta  gy  -f  Delta  a^ 

-  ■  *  ■  Delta  z"  =  Delta  gz  -f-  Delta  azr 

where  ■’■■■■■  •  ■'  ;  *'  ;  .  ■ 

Delta  x  -  Sx  x,  Delta  y  =  Sy  ~  y,  Delta  §-:■<•  z  (2.5) 
the' dfrors  of  the  inertial  systen,  J  , < 

Delta  gjj.  :  g^  Delta  gy  s  gy  ~  Syr.  Delta  gz  =  gz  feZr  '(2. 


the  errors  vhich  characterize  inaccurate  compensation  of  accelerations 
clue  to  gravitational  forces .  ,.V  '  -/  '  ■  ■  -  •  ’  ; 

If  precise  coupensation  of  accelerations  from  gravitational  s 
forces  is  accomplished  in  an  inertial  systen  (Delta  g_  »- Delta  gy  * 
Delta  gz  a  0)  and  there  is  no  drifting  of  the  gyroscopes  (Delta  * 
Delta  Oyg  *  Delta  azg  *•  0),  then  it  follows  from  equations  (2.4)  that 
the  errors  of  the  systen  will  be  determined  by  inaccuracy’  in  setting 
the  initial  conditions. 

Since  automatic  compensation  of  accelerations  fron  gravitational 
forces  cannot  be  precise  and,  in  addition,  there  is  drifting  of  the 
gyroscopes,  different  errors  do  ipecurr. In  an  inertial  systen,  and  an¬ 
alysis  of  which  is  given  below. 

3 .  Errors  in  a  Systen  of  Automatic  Compensation  for  Accelera¬ 
tions  Due  to  Gravitational  Forces.  The  errors  in  a  systen  of  this 
type  occur  because  the  compensation  signals  g  ,  g,  g ;  c,  genex-ated 
in  the  computer  are  not  equal  to  the  •components  o:r  acceleration  gx, 

Sy,  g_i  These  errors  are  fed  into  an  inertial  systen  and,  aft'erppass- 
ihg  through  a  closed  cycle,  are  again  fed  into  the  conputer.  If  the 
inertial  systen  is  stable,  the  errors  of  adtomtic  compensation  will 
be  damped.  In  the  opposite  case  they  will  grow,  which  will  give  rise 
to  an  inaccurate  determination  of  position. 

In  order  to  generate  the  compensation  signals  it  is  essential 
to  know  the  analytical  expressions  for  the  components  of  acceleration 
gx,  gy,  g2  which  are  caused  by  gravitational  forces  of  celestial  bodies 
acting  on  the  flying  apparatus,  if  tn,  Xj_,  yii  z.  (i  »  1,  2,  3,  ...,  a) 
the  masses  and  coordinates  of  the  centers  of  masses  pf  celestial  bodies 


acting  on  the  flying  apparatus,  if  tn,  Xj_,  yii  z,  (i  »  1,  2,  3,  ... 
the  masses  and  coordinates  of  the  centers  of  masses  of  ■celestial  bod 
that  are  accelerating  the  flying  apparatus;  the  components  of  the  ac 
celerations  will  be 

n  v ~ v  .  .  ( ■■■  r  .  ■*  -■  ■■■■■■',  i  ■■  ■ 


y.-Xtx 


y-y* 


Sigma  fn 
i=l  :  ‘ 


,  1 


(3.1) 


where 


Vi'  :  %  =  v  (x-Xi)2  -  (y^2 >-■  (z-z*)2  .1  ;  a  (3.2) 

and  f  is  the  gravitational  constant,  which  is  equal  to  6.67  x  10”®  cgfi. 


The  computer  should  generate  compensation  signals  gxcp  gyC,?  g2n 
in  accordance  with  the  formulas  (3*1)  and  (3.2)  in  which  it  is  neces*” 
sary  to  substitute  the  measured  inertial  ^coordinate  systen  Sx,  S  ,  S 
in  place  of  the  actual  coordinates  of  the  flyirig  apparatus  x,  y,  z. 
Consequently,  the  conputer  should  calculate  the  compensation  signals 


in  accordance  with  the  formulas 

n  S  -  x 
z  Sigma  fu.*i 


-  «<n  ®y  ~  yi 
^1X3— 


(3  -3<5 

S  -z 


Sigma  f-a 


nr? 


-  5  - 


h  --  '/'(Ex  -  nf  +  i-  <sz  -  h?  (3-4) 

In  order  to  make  the  calculations  it  is  also  necessary  to  ‘ 
know  the  masses  a.  and  the  coordinates  y^/  z^  of  the  celestial 
bodies;  which  can^be  obtained  with  great  accuracy. 

If  the  differences  Delta  g**  Delta  Delta  g„  ere  formed  by 
substituting  the  expressions  (3-1)  and  (3*3V  into  (2.6)  and  they  are 
'’decomposed  into  aeries  by  the  increments  Delta  x;  Delta  y;  Delta  z 
then/  limiting  ourselves  to  the  linear  terns ,  we  find 


Delta  g  «  Sigma 
•  x 


&.3feyj 


\  «i 


21  x-x,  (3*5) 

Delta  x-3 -g-g-  ((y-^) 


n  fin.j.4 
igma  a  W1 
i=l 

Delta  y  -f  (z-z^)  Delta  z]j* 

n  flu  /•{■  /y-y  ^  2  1  y-y, 

iff*  -~Mf  1-34— i  j  \  »8lto;y  -  2  i 

1=1  R3  il  V  R,  /  T  R  2 


Delta  g  =  Sigma 

y 


f 


[(x-xp 


Delta  x  -f  (z-zi)  Delta  z] 

n  ^i  *  V  1 

Delta  g_  #  -  Sigma  Cfl-3  ( - -  *■  ] 

■  *  ^  XRi  / 

Delta  x  4*  (y-y^  Delta  y]p 

.Substituting  (3 * 5 )  into  (2.4);  we  obtain 
■  (p2  -jh  t^2)  Delta  x  -*A  Delta  y  —  B  Delta  z  =  Delta 
—  A  Delta  x  -f  (p2  -f-  i^2 )  Delta  y  —  C  Delta  z  =  Delta  a^ 
•/)— •  B  Delta  x  —  ^.fielta^y'^-^Si2'  -^  ‘h|,fr)-Deltarz  s- Delta  azr 
where  - 

2  n  fm  i  (y-y  ,)2 

Oy.  “  Sigma  ....... . ;  il-3‘ » 

V  i=l  V  l  IL2 


z-z. 

Delta  z  —  3 — _ — [(x-x. ) 
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(3-6) 


where  • ;  .  . 

o  n  fm.  i  ,x-x, )  ] 

^  1  si«“  ft  f  3  O"  )j  ' 


(3-7) 
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Consequently,  the  errors  in  the  inertial  systen  shidh  arise 
out  of  inaccurate  automatic  coripcnsatibn  of  accelerations  due  to 
gravitational  forces  are  determined  by  a  systen  of ilinear  differential 
equations  with  variable  coefficients . ...  ,  ’  g 

Let  us  exaxiine  the  behavior  of  the  systen'  (3-6),  under1  the  con¬ 
dition  that  there  is"  no  drifting  of  gyroscopes.  The  behavior  of  the 
inertial  systen  in  respect  to  the  errors  Delta x,  Delta  y.  Delta  z 
is  detornined  by  the  type  of  solution  of  the  equations  (3.6)./rIt  i  s 
obvious  that  if  the  systen  is  stable,  the  automatic  coqpensat ion  er¬ 
rors  will  vanish,  m  case  the  systen  is  unstable,  the  errors  will  . 

grow.  .  .  .  . 

Let  us  study  the  stability  Of  the  System  of  automatic  conpen- 
sation  for  accelerations.  For  this  purpose  we  shall  examine  the 
partial  cpse,  assuaing  that  .the  coefficients  of  ;the  systen  (3-6) 
change  slowly  in  comparison  with  the  tine  required  for  processes  to 
f -take  place  in  the  autonatic  compensation  systen.  Making  use  of  this 
assumption  and  considering  the  coefficients  of  the  systen  (3-6)  to 
be  constants  for  the  sake  of  simplicity,  we  can  judge  its  stability 
on  the  basis  of  the  properties,  of  the  characteristic  determinant 


v  r 


Delta  s. 

v-;  a.  ;a;  v-*:  ..tx. 


•'  •  *  -J  ; 

-*n  B 


•  p  I  o'  ■'.*  Ofi; 
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••  ;o o  :jt:  If  this  determinant  is  expanded  .and  -set  equal  to  zero ,  we . ,  ;  >  ■). 
shall  obtain  thO  characteristic  equation  of  the  systen  (3-6) 

>  ...•••  •  '  ''  '“'hi  'c  S  ?V'i  •  -T:  (3*10) 
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• p  -mv  +  v  + v> *■* 

— C2)p2  4  \^r%E;2  “•*  A^s  — B2!^2  —C2^2  —  2ABC  8  0 


, . It  follows  fror.i  this  equation  that  the  systen  of  automatic 
condensation  of  accelerations  due  to  gravitational  forces  is  unstable  , 
as  equation  (3.10  lacks  terns  with  odd  powers  of  the  operator  p  =  d/dt. 

Lt  us  transform  equation  (3*10)/  for  which  we  shall  evaluate 
its  coefficients  by  making  use  of  formulas  (3*7)  and  (3.8).  At  the 
Sdne  tine,  we  shall  assume  for  the  sake  of  simplicity  that  the  flying 
apparatus  is  in  the  gravitational  field  of  one  celestial  body,  for  ex¬ 
ample,  the  sun.  In  case  it  is  essential  to  take  into  account  the 
gravitational  forces  of  several  celestial  bodies,  one  nay  introduce 
the  concept  of  an  equivalent  celestial  body,  for  which  g/R  n  Si§na 

WrA 

We  have  (3.11) 
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Substituting  (3.11)  into  (3.IO)  and  making  the  transfornatiohs, 
we  obtain 


where 


,  •  0  P<2  2  '  •  2v  '  •  :v:-  (3-12) 

V  ,  Cdega^r  (P  2  Onega  ;  -  0  ; 


Onega2  =  g/R 


.  f 

>•  It  follows  from  this  that  the  characteristic  equation  (3*12) 

has  two  double  imaginary  roots  and  two  real  roots,  one  of  which  is 
positive.  Consequently, .. the  general  notion  of  the  automatic  compensa¬ 
tion  systen  consists  of  three  notions :■  harmonic  oscillations  with  a 

period  of  .  - 

=  2  Pi  /  Omega  =  2  Onega  TB/g;  oscillations  with  the 

sane  period  and  an  amplitude  which  grows  proportionally  with  time;  and 
finally,  and  aperiodic,  exponentially  growing  motion  with  a  time  con- 

!  :  '  ?  =  1/ V2  Cdega  =  T,/2  V2  Pi  , 

...  ■  ,  '  •  ■  -  ‘  !  •’  '  •;  j  ;■  , 
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•'  For  a  more  detailed  examination  of -/the  character  of  the  notion 
of  the  system,  let  us  examine  the  flight  of  a  flying  apparatus  in  the 
gravitational  field  of  the  sun  in  the  direction  of  the  y^  axis,  In 
this  case  R  s  y,  x  =  z  s  0  ;  gx  =  gz  =  0,  gJ’  *  g  and  the  equations  ■ 

(3.6)  have  the  form  ■'  . ::  :i  v' 

^ 1 . ^  . ;  (3.13) 

:  i  -'  ■  (p2  4  Omega2)  Delta  x  -  Delta  a,:r  -  i  '• 

(p2  —  -2  Cmega2)  Delta  y  =  Delta  a^.  • 

.  (p2  Qroega2)  Delta  z  s  Delta  azr 

It  follows  from  this  that  in  directions 1  perpendicular  to  the 
resultant;  vector  of  the  gravitational  forces,  the  errors  Delta  x  and 
Delta  "2  change  haitionicaliy  with  the  amplitude,  which  depends  on  the 
initial  conditions  and  the  period  of  oscillation  mentioned  above.  In 
the  direction  of  the  vector  of  the  gravitational  forces,  the  error 
Delta  y  grows  exponentially.  If  there  is  a  connection  between  the ; 
channels,  as  shown  above,  a  third  notion  appears  with  oscillations  of 
increasing  amplitude.  ■;  ‘  !  -y'-  .. 

We  note  an  analogy  between  the  character  of  the  change  'Ih  the 
errors  of  the'  inertial!  system  of  navigation  and  the  notion  of  a  - 
satellite.  The  period  of  change  in  the  errors  (oscillation  period) 
of  the  inertial  system  T  \  v  :  ',:j  i  .  '  '''1'  ''  ''  "■ 

■'*_  ■■■  T  “  2  Pi  v  R/g  is  equal  to  the  periods  of  a  : 

satellite  flying  in  a  circular  orbit  with  radius  R  around  a  celestial 
body  with  the7  first  cosmic  velocity  .  ‘  '  :v'.  ' j  r  ;- ' 

-■  V^;=  R  Onega  “  V  gR;  (3 *  I1*- ) 

When  navigating  at  the  surface  of  the  earth  the  period  of 
change  in  the  errors  Delta  x  and  Delta  z  is  equal  to  =  84.4  minutes 
and  is  called1  the  Schuler  period.  '■  ' 

1  The  error  in  the  automatic  Compensation  of  accelerations  duo 
to  gravitational  forces  Delta  y  in  the  direction  of  the’  resultant  vec¬ 
tor  of  the  gravitational :  forces'  changes  exponentially  with  the  constant 

of  title  •,  _ , -  i _ _  :'--J 

T  =■  4- *\lR/g,  which  is  l/2  V  2  Pi  of  the  period  of  the' in-, 
ertial'  hyiten.  ^Wheh  the'  initial  'cohditio&'s  -cat  nonzero)  the  error  ,'V 
grows ‘ continually,  which  characterizes  the'  instability  of  the  automatic 
compensation  system.  ’  Consequently,  switching  on  the  automatic  compen¬ 
sation  channel' of  an  inertial  system  in  which  the  axis  of  the  accelero¬ 
meter  coihcides:  with 'the  direction -of  the  vector  of  the  gravitational 
forces  takes  the  system  absolutely  unstable. f.  i  1  i- :  ; 

;  The  analogy  betiieeii  the  inertial  system  and  the  notion  of  a 
satellite  can  be  continued.  We  note  in  particular  that  the  character 
of  the  change  in  the  error  Delta  y  corresponds  with  the  change  in  the 
radius  vector  R  of  a  satellite  which  is  moving  away  from  a  celestial 


body  (for  example,  fron  the  earth)  with  the  second  cosuic  velocity 

_ _  ,  (3.15) 

V2  =  v  2  Onega  R  =  T2gR 

Indeed,  only  when  the  flying  apparatus  (satellite)  is  moving 
at  the  second  cosmic  velocity  does  it  overcome  the  attraction  of  the 
celestial  body  and  move  away  from  it  freely.  In  the  inertial  system., 
the  signal  in  the  automatic  compensation  channel  which  measures  the 
coordinate  Delta  y  fully  compensates  for  the  acceleration  of  gravity; 
therefore  the  inertial  system  is  operating  under  conditions  analogous 
to  the  physical  absence  of  gravitational  force. 

The  analogy  between  the  motions  of  an  inertial  system  and  a 
satellite  can  be  used  for  simulating  the  motion  of  a  satellite. 

In  order  to  evaluate  the  effect  of  the  initial  error  Delta  y q 
and  the  altitude  of  flight  above  the  earth  on  the  magnitude  of  the 
error  of  the  inertial  system  Delta  y,  determined  by  means  of  the 
second  equation  of  (3*13)  with  a  zero  right-hand  side,  the  correspond¬ 
ing  graphs  of  Delta  y/Delta  yQ  =  f(t)  are  presented  in  Figure  5*  Die 
curves  1,  2,  and  3  corresponduto  flight  altitudes  of  0,  1,000,  and 
2,000  kilometers.  The  error  accumulates  more  slowly  with  increased 
altitude,  since  the  effect  of  gravitational  forces  decreases. 

The  properties  of  an  inertial  system  which  were  studied  here  . 
correspond  not  only  to  the  special  case  of  flight  along  a  straight  1 
line,  but  also  to  arbitrary  flight  if  the  axis  of  sensitivity  of  one 
of  the  accelerometers  is  held  in  the  direction  of  the  vector  of  the 
gravitational  forces.  Under  terrestrial  conditions,  this  means  di¬ 
recting  the  axis  of  sensitivity  Of  one  of  the  accelerometers  along  the 
vertical  direction  and  the  axes  of  sensitivity  of  the  other  two  ac¬ 
celerometers  in  the  horizontal  plane.  The  errors  in  the  channel  which 
measures  the  altitude  of  flight  will  grow  exponentially  with  zero  in- 
itial  conditions,  while  the  errors  in  the  channels  which  measure  the 
coordinates  of  position  will  grow  harmonically  with  the  M.  Schuler 
period. 

When  a  flying  apparatus  moves  along  an  arbitrary  path,  all  the 
channels  of  a  system  of  automatic  compensation  for  accelerations  due 
to  gravitational  forces  are  unstable,  as  shown  by  equations  (3*6). 

When  the  initial  conditions  are  nonzero,  the  errors  in  all  three  chan¬ 
nels  have  terms  which  grow  with  time  in  addition  to  the  oscillatory 
terms. 

The  instability  of  inertial  systems  with  automatic  compensation 
for  accelerations  due  to  gravitational  forces  cannot  constitute  a 
criterion  of  their  inapplicability.  Errors  will  not  exist  under  zero 
initial  conditions,  and  in  case  initial  errors  do  exist,  their  growth 
is  slow.  In  this  case,  the  time  of  growth  is  comparable  with  the 
period  of  the  system.  Moreover,  the  system  can  be  made  stable  by  in¬ 
troducing  additional  links  or  by  feeding  additional  external  informa¬ 
tion  into  the  system  [2],  for  example,  from  astronomical  navigation 
systems.  -  .■  -  ... 
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has 


One  can'  f ohnulate -the  following  theorem  on  the  basis  .of  -What 

been  set  forth  here.  * 

Theorem  1.  Automatic  conpensation  for  accelerations  due  to  ; 
gravitational  forces  in  all  three  channels  of  an  iiiertial  system  will 
lead  to  the  characteristic  equations  .. 

(3-ib; 

' q' ;  1  (p2  ttaega2)2(p2 .?  Onega2)  *  0  ■■■,<  ,<?  ... 

for  the  errors  in  each  of  the  channels. 

Theorem  2.  The  period  of  oscillation  of  the  errors  of  ^.n  .in¬ 
ertial  system  is  equal  to  the  period  of  a  satellite  moving  in  a  circu¬ 
lar-  orbit  with  radius  R  around  a  celestial  body  at  the  first  cosmic 

-velocity  V,  •-■‘V.'gK.  ...  ;j.;,  -  ,,, ; 'rl'i  i 

.>v  c, Theorem  3.  'gte  constant  of  the  tine  of  growth  of  ttie  errors 
of  an  inertial  system  is  equal  to  the  constant  qf  the  tine  of  motion 
of  a  satellite  which  is  moving  away  from  a  celestial  body  with  the^ 
second 'cosmic::  velocity  :'.h- so w’’  .-i?-"' 

'■'■"■•f-  ;in  srcri  f.  off  'io  a .Vp:  -*  y  .2sR«  ...:n  t  f  ■  .. 

OA  method  of  "freezing";  the  coefficients  is  used  in  studying  the 
equations'  (3.6 )  and  as  a  result  we  obtain  a  system  With  constant  Coef¬ 
ficients.  This  assumption  can  be  considered  adnissable  fpr small  in¬ 
tervals  of  time  .comparable  with  the  period  of  the  oscillations  of  an 
-inertial  system -and  the  conclusions  given  above  are  .justified.  For 
; large  intervals  of  time  it  is  essential  to  examine  equations  with 
variable  coefficients .  ,-In  this  case  the  quantities  Tx  and  T2  become 
variables, --'thus  ’the  period  of  oscillations  and  "the  constant  of  time 

lose  their  usual  sense;.,.;  : : ./:+  '  '  . 

.u4v  rJnstrumental  Errors .-  Instrumental  errors  in  ati  inertial 
system  consist  in  drifting  of  the  gyropiatf ofm  relative  to  ~  the  iner¬ 
tial  coordinate  system,  errors  in  the  accelerometers,  the  integrators, 
the  computer,  et  cetera.  ^Rrifttag  ;of '  ihe  gj^cplatfo^  of 

errors  :.of  .  two  types i  ;  the  errors  ;  linked,  with  rotation  of  the  xyz 
coordinate  system  relative  to'  the  inertial  XpygZQ  coordinate  system 
and  the  errors  generated  in  a  closed  compensation  system.  :; 

An  idea  .of  the  effect  of  instrumental  errors  on  the  accuracy 
•  of  an  ihertial  system  Can  be'  obtained  from  the  following  considera- 

;<8  o.i  .  For  Delta,  -a.,;,  .>  Delta  bx^>  .  .  . ,  Delta  bx2,  ;.  .  which  are 

the  errors  of  -  the-  accelerometer,  the  first  and  second-  integrators 
respectively:,  equation  (2 .3 )  which  determines  the  calculatGdcoordinates 
acquires  the  form  '  >  'rA.  .  f .  ,i  ■  ! .. 1 ' (4.1) 


" .  .  .  (  t  T  + 


I 


l(x"  -f  Delta'  g£  -j-:  Delta  axr  'Hh  Delta  a^)  --l-f  Delta  bxl 

L  .  -  ■  •  :  7!J  :r K  ::  \-  w*  -sn  }*  "V'.:* 

k_  _ _  .  1  ..  „  ’ 


+  v) 


2x  -fr  Delta  bx2  -f  xQ  =  S, 
p 


Continuation  of  (4.l) 

It  i 

[(y"  -f  Delta  Sy  4  Delta  nyr  4  Delta  ayG)  "jr  Delta  byl  -f  yQ’ j 

.  4  Delta  b  ^  4  yn  -  S„  , 

■"~p —  ■  y2  °  -  y  •  ;  -■ ;  - 

r  ?  kiz  ] 

j(z"  4  Delta  gz  4  Delta  azr  -f  Delta  aZQ)  4  Delta  bgl  4  zQ'  ] 

k_  >  ,  >  ....  •  *  V:  I"- .  ’  r'  v‘  ;  '  ;. . 

£z  .4  Delta  bz2  rf-  zQ  =  Sz  ;  ;  ■ 

P  It  follows  froa  this  that  instrunental  errors  in  the  elements 
of  an  inertial  systen  have  a  direct  effect  on  the  accuracy  with  which 
the  coordinates  of  position  are  determined.  Instrunental  errors  will 
reach  large  magnitudes  when  an  inertial  systen  is  operated  over  pro¬ 
tracted  periods  of  tine.  -  !  ■  ' 

The  principal  nethods  for  reducing  instrunental  errors  are 
essentially  reducing  the  errors  of  the  elenents  of  the  inertial  systen 
,  and  correcting  the  systen  by  introducing  external  navigational  infor¬ 
mation.  In  the  latter  case  the  corrections  can  be  made  discretely, 
not  continually.  '  s:  .  .  /  _ 

:  .  In  concluding,  we  note  that  an  Inertial  systen  With  three  chan¬ 
nels  of automatic  conpensationffor  accelerations  due  to  gravitational 
forces  is  unstable  when  a  flying  apparatus  is  moving  arbitrarily  and 
that  the  errors  in  this  systen  will  grow  When  the  initial  conditions 
are  nonzero .  In  order  to  eliminate  the  errors  it  is  necessary  to  c  . 
correct  the  systen  continuously  or  discretely  by  feeding  external 
.  navigational  information  into  it.  In  this  case  the  inertial  systen 
,  fulfills  the  role  of  a  memory  device  in  addition  to  its  basic  role. 

■  The  analogy  between  the  changes  in  the  errors  of  an  inertial 
system  with  three  automatic  compensation  channels  and  the  notion  of 
•  a  satellite  can  be  used  in  simulating  the  notion  of  a  satellite. 
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